%% 循环结构

% 一般格式
% for 循环变量=矩阵表达式
%     循环体语句块
% end

% 执行过程是依次将矩阵的各列元素赋给循环变量，然后执行循
% 环体语句，直到各列元素处理完毕循环结束

% 利用不定积分求解圆周率π
% 用不定积分求解四分之一圆的面积
n = input('分割成多少个梯形：n=');
h = 1/n;
x = 0:1/n:1;
y = sqrt(1-x.*x); % 不定积分的被积函数
s = 0; % 梯形面积累加和
for i = 1:n
    s = s+ (y(i)+y(i+1))*h / 2;
end

pi = s*4;
disp(pi);

